Velocity selection problem in the presence of the triple junction 
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Melting of a bicrystal along the grain boundary is discussed. A triple junction plays a crucial 
role in the velocity selection problem in this case. In some range of the parameters an entirely 
analytical solution of this problem is given. This allows to present a transparent picture of the 
structure of the selection theory. We also discuss the selection problem in the case of the growth of 
a "eutectoid dendrite" where a triple junction is present because three phases are involved in the 
eutectoid reaction. 
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During the last decades, our understanding of pattern 
formation in various nonlinear dissipativc systems has 
made a remarkable progress. Building on these founda- 
tions, it has now become possible to develop a descrip- 
tion of a large class of patterns that are found in diffu- 
sional growth. In the two-dimensional dendritic growth 
the needle crystal is assumed to be close to the parabolic 
Ivantsov solution If anisotropic capillary effects are 
included, then a single dynamically stable solution is 
found for any external growth condition. In this classical 
problem of dendritic growth, velocity selection is con- 
trolled by tiny singular effects of the anisotropy of the 
surface energy (see, for example, Q, [I])- In the case of 
isotropic surface energy, the dendritic solution does not 
exist and instead the so-called doublon structure is the 
solution of the problem 0, 0| . The dendritic-like struc- 
tures which occur during melting have also been observed 
experimentally (see [(§ and references therein). 

Moreover, for many applications, such as grow of 
dislocation- free silicon crystals [3], the investigation of 
the influence of the defects on the melt-crystal interface 
behavior is of great importance. The melting along de- 
fects such as grain boundaries in inhomogencous mate- 
rials is typically favorable, since they serve as predeces- 
sor to melting. In this paper we therefore discuss the 
melting process along the grain boundary in the crys- 
tal (see Fig. 1). The presence of a triple junction at 
the tip of the melting zone leads to an entirely different 
selection mechanism, since the triple junction produces 
a very strong perturbation of the liquid-solid interfaces 
and weak anisotropy effects can be neglected. This idea 
was expressed in [8|. We also discuss importance of a 
triple junction for the selection problem of the growth of 
a "eutectoid dendrite" (see Fig. 2). 

Melting along a grain boundary. We consider the two 
dimensional melting problem of a pure bicrystal with a 
straight grain boundary as shown in Fig. [T] Overheating 
of the crystal leads to a moving melt zone. The contact 
angle <f> at the triple junction is related to the surface 
energies by Young's law: 27 cost/) = 7&, where 7 is the 
surface energy on the liquid-solid interface and jb is the 
surface energy of the grain boundary. We introduce the 




FIG. 1: Propagation of the melt zone along the grain bound- 
ary (x — 0). The melting front propagates with velocity v 
along the y-direction. S corresponds to two solid grains and 
L corresponds to the liquid phase. The triple junction is at 
the origin of the coordinate system. 



dimensionless temperature field u = c p (T 00 — T) / L, where 
L is the latent heat, c p is heat capacity, and is the 
temperature in the crystal far away from the interface. 

The temperature field u obeys the following heat dif- 
fusion equation and boundary conditions: 



DV 2 u = du/dt, 
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where L and S refer to the liquid and solid, respectively, 
c?o = jTmc p /L 2 is the capillarity length, Tm is the melt- 
ing temperature, and k is the curvature of the liquid-solid 
interface, which assumed to be positive for convex inter- 
faces (as in Fig. 1). We introduce the dimensionless 
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FIG. 2: Schematic picture of the interface structure in the 
eutectoid reaction. The structure propagates with a constant 
velocity v in the y direction. The triple junction is located at 
the origin of the coordinate system. 



overheating A = c p (T OQ —Tm)/L, and the Peclet number 
p = vR/2D with D being the thermal diffusion constant 
(assumed to be the same in both phases), R is the radius 
of curvature of asymptotically fitted Ivantsov parabola 
and v is the velocity of the steady-state motion. 

The physics underlying Eqs. ([TJS]) is quite simple. A 
moving melting front absorbs the latent heat. Require- 
ment of the heat conservation at the interface gives us 
Eq. ([2]) (n is the normal to the front pointing into the 
solid phase; v n is the normal velocity of the front). The 
local thermodynamical equilibrium in the interfacial re- 
gion implies Eq. Q. This is the Gibbs-Tomson relation, 
which gives the equilibrium value of the temperature at 
the interface taking into account curvature corrections. 

An equivalent formulation of the problem which is 
more convenient for our purposes is obtained by elim- 
inating the thermal field. This can be done by using 
the standard Green's function techniques to obtain an 
equation in closed form for the shape of the solid-liquid 
interface: 
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where i] = [(x — x') 2 + (y — y') 2 ] % , and Kq is the modified 
Bessel function of third kind in zeroth order; all lengths 
are measured in units of R. The relation between the 
overheating A and the Peclct number p is given by the 
Ivantsov relation: A = v /p7rexp (p) erfc(y^). 

Eq. ([4]) is a complicated nonlinear intcgro-diffcrential 
equation. We should find a smooth solution of this equa- 
tion which has a proper angle at the tip and which is 
close to the Ivantsov parabola (y = —x 2 /2) in the tail 
region. The classical dendritic growth problem does not 



have a solution with isotropic surface tension and with a 
smooth tip which corresponds to the angle <f> = tt/2 @,H|. 
This statement can be expressed in the following form. 
For any given positive values of the Peclet number p and 
the so-called stability parameter tr = da/pR the symmet- 
ric solution which is close to the Ivantsov parabola in the 
tail region has an angle at the tip = f(cr,p) < tt/2. The 
limit cr = and <fi = tt/2 is a singular limit for that prob- 
lem. For example, Meiron 0] numerically calculated the 
angle <j> as a function of a for several values of the Peclet 
number with isotropic surface tension and found that the 
angle (j> remains smaller than tt/2 for any positive a. 

Our problem differs from the classical dendritic prob- 
lem precisely in that we need to satisfy the condition 
that the angle <fi at the tip should be smaller than tt/2. 
This means that we can select the stability parameter 
a = a*(<p,p) as a function of and p. In order to illu- 
minate this selection mechanism we consider the case of 
small opening angles i.e., f)« 1 and small overheating 
(p <C 1) which allows us to drastically simplify the prob- 
lem and finally obtain the selection in some limit entirely 
analytically. 

For small opening angles </> the following rcscaling of 
coordinates is useful in order to eliminate small parame- 
ters from the interface shape: 

x^x/cj), y^y/4> 2 . (5) 

After this rescaling, y = —\x\ near the tip and in 
the tail region the asymptotic of the Ivantsov parabola 
(y = —x 2 /2) remains unchanged. In this small angle ap- 
proximation the function n depends only on y variable, 
r\ ss \y — y'\(j> , because x <C y. Consequently the whole 
integral kernel in Eq. (|4]), under the assumption of small 
Peclet numbers, is function only of y. Due to the sym- 
metry of the interface we consider it only for x > 0. It 
is more convenient to treat the shape as a function x{y) 
and to make change of variables, 

y -> -y- 

The curvature of the interface can be written as n ss 
—d 2 x/dy 2 . Performing all these step we reduce the orig- 
inal nonlinear Eq. (j4]) to the following linear equation: 

d 2 x 2 I p f°° , ,dx(y') -'(,'-«)„ , p , , ,. 

(6) 

where fid = do(j) 3 /RA = a^^fpfrr. As follows from 
Eq. ([6|), the eigenvalue [id is a function only of one pa- 
rameter, p/4> 2 . We solve Eq. © numerically using stan- 
dart linear algebra routines. The results are presented in 
Fig.H 

For small values of p/(f> 2 <C 1, the eigenvalue fid is pro- 
portional to the prefactor in front of the integral term in 
Eq. © i.e., fid = a^J (p/4> 2 ). We found numerically that 
a sa 2.2, which leads to the following scaling relations for 
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FIG. 3: Eigenvalue /id as a function of p/4> 2 ■ Solid line repre- 
sents the numerical solution of Eq. (|6]) , dash line corresponds 
to the asymptotic behavior fid = 2.2^p/ <f> 2 . The inset shows 
an enlarged region of small values of p/4> 2 , the dashed line 
is 2.2^/p/ '(f> 2 , crosses correspond to the data extracted from 
Ref.0. 




FIG. 4: The schematic phase diagram of the eutectoid reac- 
tion. To is the temperature of the system, Te is the eutectoid 
temperature. The regions 7, a and j3 correspond to the one 
phase equilibrium state. 7 + a and 7 + /3 are regions of two- 
phase equilibrium. C' a , Cp and Ce are the equilibrium con- 
centrations at the eutectoid temperature. For this diagram 
the eutectoid composition Ce = (C a + C[j)/2; the lines of 
the 7 — a phase equilibrium are parallel, the same for lines of 
7-/3 equilibrium; the equilibrium concentrations at the a — f3 
interface, C a and Cp do not depend on the temperature. 



the selected tip radius and velocity: 
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We compare our asymptotic predictions fid = 2.2y / p/</> 2 , 
with data extracted from Fig. lc of the Ref. @. Despite 
the fact that the results of Ref. [s| were obtained by the 
solution of the full Eq. [4J and our results were obtained 
under the assumption of small opening angles, they show 
a relatively good quantitative agreement up to angles <p ~ 
7r/6 (see the inset in Fig. [3]). 

In the case of large value of p/(j> 2 1 further sim- 
plifications of Eq. ([6]) are possible. The combination of 
the exponential term with the asymptotic of the Bcsscl 
function for the large arguments leads to the integral ker- 
nel 1/ \/p(t>~ 2 (y — y') for the argument y' < y, while, for 
y > y' the integral kernel is exponentially small and can 
be neglected. Thus, Eq. ^ finally reads 
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This equation for the unknown function f(y) = dx/dy 
can be solved by means of the Laplace-transformation 



F i s ) = f{y)exp(sy)dy. 
Jo 

Using the boundary condition, /(0) = 1, we get 



F(s) 
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3/2 _ 



(9) 



The function F(s) has a pole at s = (2/irfi 2 ,) 1 / 3 . There- 
fore the physically relevant solution exists only if the pa- 
rameter fid = 7r/2, which means that the pole and zero of 



the function F(s) coincide. To our best knowledge this is 
the first example where the problem of velocity selection 
in dendritic growth is solved entirely analytically. This 
exact analytical treatment demonstrates in a very trans- 
parent way why the physically relevant solution of Eq. @ 
exists only for some specific value of the parameter \id- 
The inverse Laplace-transform of Eq. © finally gives the 
closed expression for the derivative of the interface pro- 
file: 



dx 
dy 




x crfc 



We note that the same analytical treatment is possi- 
ble also in the case where interface kinetics effects arc 
included. The interface kinetics requires an additional 
term, —v n //3, in the left-hand-side of Eq. (j4]) and, corre- 
spondingly, the term —\ipdxjdy in Eqs. [8j) (J5 is the 
kinetic coefficient and = v<fi/((3A)). In this case the 
selection condition reads fid = tt(1 — M/3) 2 /2- 

Eutectoid "dendrite " 

Another illustrative example of selection by the triple 
junction is the eutectoid reaction. We consider here a 
rather simple model of the eutectoid system, described by 
the phase diagram in Fig. [51 in order to obtain the equa- 
tion in the spirit of Eq. ([4]). We restrict the consideration 
to the case of an isolated eutectoid dendrite (Fig. 2), i.e., 
we consider the situation when the period of the lamel- 
lar structure is much larger than the characteristic tip 
scale of the "dendrite". We call the appearing structure, 
shown in Fig. [51 "dendrite" because, as we show later, it 
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FIG. 5: Eigenvalue a as a function of the opening angle. The 
solid line (<f) > tt/2) corresponds to the growth of a and (5 
phases from the 7 phase, To < Te- The dotted line (<j) < 7r/2) 
to the inverse process, To > Te- 



has the parabolic asymptotics. The concentration fields 
have jumps on the phase boundaries in cutcctoid reac- 
tions. As a consequence the difference of concentration 
Ce — C a (source term), and Ce — Cp = {Ce — C a ) (sink 
term) appear as the amplitudes in front of the integral 
terms. 

Moreover, we also assume that the diffusion constants 
are equal in all three solid phases. This assumption cru- 
cially simplifies the model. The concentration C 7 far 
away from the interface is chosen to be Ce leading to the 
symmetry of the shape as presented in Fig. 2. Finally, 
this system can be described by the integro-differential 
equation as in the case of classical dendritic growth. Due 
to the assumed symmetry of the phase diagram and of 
the shape of the interfaces, we consider only positive val- 
ues of x, i.e., write down the equation only for the 7 — ft 
interface: 



* d p 

A K = 

R TX 
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dx'e- p (y-y">K (pr)(x,x')) 

O 

dx'e- p(v - v,) K {pr](x,x')). 



(10) 



Here do is the capillary length [n|, k is the curva- 
ture of the interface which is positive in Fig. 2, and 
A = {Ce — C?Tp)/(Cp — Ce) is the supersaturation. This 
equation differs from Eq. (j4|) only by the sign in front 
of the first integral, which represents the source at the 
a — 7 interface, while the second integral corresponds to 
the sink at the (3 — 7 boundary. In other words, this 
simple modification of the classical equation of dendritic 
growth already contains the key ingredient of the eutec- 
toid reaction. 

In the most relevant case of small Pcclet numbers this 



equation can be further simplified. We stress that the 
usual relation between the Peclet number and the super- 
saturation A w ^firp , is no longer valid here. In the 
small Peclet number limit the flux in the a and (3 phases 
is much larger than in the 7 phase. Thus, in the asymp- 
totics, one can write DA/x = —vdx/dy. This leads to 
the parabolic asymptotics of the shape, y = —x 2 /2R, and 
to the following relation between the supersaturation and 
the Peclet number, p = vR/2D: 

A = 2p. 

Thus, for small supersaturation A Eq. (|10[) reads: 



9 1 r . ( {x + x'f + {v-y') 2 



where a = do/Rp. Our purpose is to find the shape 
with parabolic asymptotics and the eigenvalue a for the 
given angle <p. We solve this nonlinear equation numer- 
ically in the spirit of Ref 0]. The resulting eigenvalue 
a" = c*(</>) as function of the opening angle is shown in 
Fig. [5] We notice, that in contrast to the problem of clas- 
sical dendritic growth, where a = c*(0,p) has a singular 
point at cj> = 7r/2, in the case of the eutectoid reaction, 
a solution exists even for (j> = tt/2. Angles <fi > 7r/2 cor- 
respond to growth of a and (3 phases from the 7 phase 
(To < Te), while angles cf> < tt/2 describe the inverse 
case (To > Te). In the latter case the structure of the 
interfaces is the same as presented in Fig. [T] with the 
replacement of L — > 7 and S — > a, f3. The related model 
of melting of eutectic structures was discussed in [11 . 

Since the eigenvalue a* is found, the velocity of the 
process can be written asi> = 2a* Dp 2 / do. Consequently, 



the velocity of eutectoid growth is proportional to A 2 , 
while for the classical dendritic growth it scales as A 4 . 

In conclusion, the selection mechanism due to the pres- 
ence of the triple junction is very different from the 
anisotropy effects which are responsible for the selection 
in the classical theory of dendritic growth. In the regime 
of small opening angles the problem of velocity selection 
in melting along the grain boundary has been solved en- 
tirely analytically. This exact analytical treatment rep- 
resents in a very transparent way the structure of the 
selection theory. We hope that results of this paper will 
stimulate some new experiments to bring new insights to 
this interesting problem. 
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